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ABSTRACT 

Spectral energy distribution (SED) in protostellar disks is determined by the disks' 
internal dissipation and reprocessing of irradiation from their host stars. Around T 
Tauri stars, most mid-infrared (MIR) radiation (in a few to a few ten fim wavelength 
range) emerge from regions around a fraction to a few AU's. This region is interesting 
because it contains both the habitable zone and the snow line. Recent observations re- 
veal that SED variations, in the MIR wavelength range. These variations are puzzling 
because they occur on time scale (a few days) which is much shorter than the dynam- 
ical (months to years) time scale at lAU to a few AU's. They are probably caused by 
shadows casted by inner onto outer disk regions. Interaction between disks and their 
misaligned magnetized host stars can lead to warped structure and periodic SED mod- 
ulations. Rapid aperiodic SED variations may also be induced by observed X-ray flares 
from T Tauri stars. These flares can significantly modulate the ionization fraction of 
the gas and the net charge carried by the grains near the surface of the inner disk. The 
newly charged grains may be accelerated by the stellar or disk magnetic field and adjust 
their distances from the midplane. Shadows casted by these grains attenuates the flux 
of stellar photons irradiated onto regions at several AU's from the central stars. We use 
this model to account for the observed rapid aperiodic SED variabilities. We suggest 
regular monitoring of SED variations will not only provide valuable information on the 
distribution of the disk aspect ratio near the habitable zone but also provide a probe on 
the interaction between the inner regions of the disk with the magnetosphere of their 
host stars. 



1. Introduction 

For mid-infrared spectrum in protosteller disks, normally we expect that the spectral variation 
is determined by the local disk dynamics, with typical timescale comparable to the orbital period. 
Spectrum around a few (< 10) ^um comes from disk irradiation at 0.1 ~ 1 AU, with a typical 
dynamical timescale as 0.1 ~ 1 year. However, a series of recent observations of T Tauri stars 
by Spitzer revealed the rapid variability of spectral energy distribution in protostellar disks with 
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timescales as short as one week . Muzerolle et al. sh owed remarkable mid-infrared variabihty of 



the transitional disk LRLL 31 ([Muzerolle et al. 2003 ). and highlighted the seesaw feature: "the 



variability observed at shorter wavelengths is reversed at longer wavelengths, with a constant pivot 
point at A ~ 8.5/im." The largest flux difference reported is 60% at short wavelengths and 30% at 
long wavelengths, and an impressive change (20%— 30%) occurred over just one week. Moreover, the 
flux variations are fairly constant as a function of wavelength within intervals 5 — 7/um and 11— 35/um, 
respectively. As a prevailing ph enomenon, variab ility of disk emission has also been reported in pre- 



main sequence Herbig Ae stars (jSitko et al. 20081 ). Structural changes in the region of the disk near 



dust sublimation zone were proposed to account for the variability. For relatively slow variabilities, 
Bary et al. (2009) reported significant variations of 10/im silicate features on monthly and yearly 
timescales. They further pointed out several scenarios that might produce such variability, including 
intervening cooler dust grains, disk shadowing and illurn ination, turbulent mixing, and disk winds. 



Most recently, an extensive study (jFlahertv et al. 201 ll ) has been carried out for the transitional 



disk LRLL 31 with multiple epoches of observations. Based on the observations, authors reasonably 
put some constraints on the physical mechanism responsible for the variability based on their 
observations: "the variable accretion and wind are unlikely to be the cause of the variability, nor 
is a companion within OAAU perturbing the disk." The most likely explanation they provided 
is "either a companion beyond ~ OAAU or a dynamic interface between the stellar magnetic field 
and the disk leading to a variable scale height and/or warping of the inner disk." 

Following these puzzling observations, several mechanisms have been proposed to explain the 
rapid mid- infrared varia bility. Flaherty and Muzerolle considered non-axisymmetric structure vari- 



ations in the inner disk (jFlahertv and Muzerolle 2010l ). such as a warp or a spiral wave, and sug- 



gested the observed variability could be explained by a warped inner disk with variable scale height. 
As far as we know, this is the only paper trying to construct a physical model to explain the rapid 
mid-infrared variability at present. Though their model can generate the pivot point in the SED 
with a consistent timescale, their the location of the pivot point is inconsistent with the obser- 
vations. Moreover, their model doesn't fit the flux change at A > 8fim. Based on the CoRoT 
light-curves morphology of 83 previously known classical T Tauri stars, Ale ncar at al. associate d 



the light curve variations with a magnetically controlled inner disk warp (jAlencar et al. 2010l ). 
They further pointed out the inner warp dynamics should result from the interaction between the 
stellar magnetic field and the inner disk region. However, they haven't built up a physical model 
to consolidate their proposition , and their explanation concludes a timescale of a few rotational 
periods, which usually extends beyond one week, thus inconsistent with the observations. Espaillat 
et al. (2011) showed that they can reproduce the observed variability using irradiated disk models 
by changing the height of the inner disk wall by ~ 20%. Their model builts a reasonable link 
between the SED variability and the variation of the inner disk structure. However, they haven't 
provided a physical explanation about how inner disk wall can change in a short timescale. Based 
on a self-consistent physical model, we show in this paper, that observed X-ray flares can account 
for the rapid mid- infrared variability in protosteller disks. 
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Stellar X-ray flares are prevailing among Young Stellar Objects (YSOsl (iFeigelson and Montmerle 19991 : 
Glassgold et al. 20051 : iFavata et al. 20051 : iGetman et al. 20081 : IStelzer et al. 20051 ) , which may come 
from material transfer from disk to the central star, or stellar magnetic reconnections. A recent 
COUP survey of the Orion Nebula Cluster (with Chandra space telescope) indicates that nearly 
all T Tauri stars emit X-ray (0.5-8 keV) photons with a nearly constant quiescent flux in the range 
of '~ 10^^ J s~^. D uring COUP's tw o week observing period, 1-2 flares are commonly observed 
around each source ijWolk et al. 20051 ). These flares are characterized by power-law energy spectra 
and up to 1-2 orders of magnitude increase in their peak X-ray luminosity and they decay expo- 
nentially on a time scale ~ 10'^~^s. Simil ar flares (some with larger amplitudes) were observed with 



the X EST survey of the T aurus region ljCudel and Telleschi 20071 ) and the Serpens star-forming 
region ljGiardino et al. 20071 ). They have be en modeled with e mission from thermal plasma along 
cylinderical loops which extend beyond fim ljReale et al. 19971 ). These X-ray flares can greatly in- 
fluence the circumstellar disks in various aspects. For example, (jWolk et al. 20051 ) estimated that 
the X-ray can ionize protostellar disks out to lAU at a rate which is an order of magnitude higher 
than that estimated for cosmic rays. In this paper, we focus on the ionization effect of the X-ray 
flares. 

We develop an analytical framework to show that the observed rapid X-ray flares from T Tauri 
stars can significantly modulate the ionization fraction of the gas and charges carried by the grains 
near the surface of the inner regions of their surrounding disks. In the proximity of their host stars, 
the newly charged grains are accelerated by the stellar and disk field and adjust their distances 
from the midplane. Consequently, the shadow casted by the grains on the surface of the inner disk 
regions attenuates in the flux of stellar photons irradiated at several AU's from the central stars. 
We use this model to account for the variabilities of spectral energy distribution, at a few micron 
range, on time scales of a few days. 

The rest of this paper is organized as follows. Our primary objective is to identify and ana- 
lyze the dominant physical process rather than to reproduce the SED evolution with the detailed 
construction of a new set of disk structure models. In Section 2, we briefly describe the static 
disk structure, mainly based on a model of DuUemond et al. (2001). We analyze X-ray ionization 
and grain charge loading process in the inner frontier of the disk in Section 3, and calculate the 
grain dynamics in Section 4. In Section 5, we formulate the spectrum variation, and discuss the 
response-time to verify the viability of our analysis. In Section 6, we present the main features 
of our model through a numerical example, with comparisons with observations. Finally, in sec- 
tion 7, we conclude the paper and discuss different scopes in parameter space to further explore 
implications of our model. 



2. Disk Model 



In their seminal paper, Chiang and Goldreich (1997) proposed a hydrostatic radiative equilib- 
rium model for passive disks surrounding T Tauri stars. The most notable idea they developed is a 
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"sandwich" model for the vertical structure of the disk, which consists of an interior layer amidst 
two surface layers. The optically thin surface layer is directly heated by the stellar irradiation, and 
the interior region is heated by the re-irradiation from the surface layer. The direct implication of 
this model is the flaring configuration of the disk. Dullemond et al. (2001) modified this model 
by adding an inner hole due to the evaporation of carbon or silicon dusts in the disk. In their 
model, the inner rim of the disk is puffed-up, because of direct exposition to the stellar flux. As a 
consequence, there is region right beyond the inner rim which is in the shadow. Garaud and Lin 
(2006) further extended the disk structure by considering the both internal dissipation and surface 
irradiation. 

In this section, we generally follow the line of Dullemond et al. (2001) with some modifications 
to incorporate internal dissipation, which dominates the energy spectrum in shadowed region. For 
sake of simplicity, we write down the main results of their works without detailed discussions. We 
consider a T Tauri star with a circumstellar disk, which is heated by both the irradiation of the 
central star and the internal viscous dissipation. Under the framework of the Minimum-Mass Solar 
Nebula Model, the surface mass density can be expressed in the form S = 'Eo{R/AU)~^^'^ , with 
Sq = IQ^kg/m?. The disk is divided into several regions: the puffed-up rim which is the inner 
edge of the disk locating in the sublimation line of silicates; the shadowed region which lacks the 
direct irradiation of the central star by the rim sheltering; the flaring part including the interior 
and surface layers. The intermediate shadowed region is heated by viscous dissipation whereas the 
outer part is mostly heated by stellar irradiation. For protostellar disks with a solar composition 
and accretion rate around solar mass host stars, the transiti onal radius between these two 



regions is nran ~ l(Md/3 x IO^^Mq^-^)"-^ (|Garaud and Lin 200fil ^ 



In most parts of the disk, small (/xm-size) grains are thought to be homogenously mixed with 
the gas, with the mass density percentage Z = Bmi^iHL = i%_ w/g take grains to be spheres with the 
average radius ranging from 0.1^?n to Ifim, and the density ~ 3 x lO^kg/m^ . The emissivity 
and correspondingly opacity possess two Gaussian peaks located at 10pm and ISfim due to silicate 
resonance. When dealing with the heat equilibrium equations, we omit the contribution of silicate 
resonance, and use e = T/T^, and k = k^T/T^ instead. 

Optical height H and temperature in the interior layer Tj at radius R can be expressed as 



7 J \2'iljiJ \ToJ \R 



7J ^ \2^J \ToJ \R 
The dimensionless constant x varies between 2 and 6 in the region discussed, w hich can be deter- 



mined with an appropriate deflnition of the optical height (jGaraud and Lin 20061 ) . Functions ipi and 



tps depend on the temperature in the interior layer and surface layer respec tively, accounting for th e 



possibility that the disk layer is not fully optically thick to its own emission (jPullemond et al. 200ll ) 



The symbols and i?* represent, respectively, the surface temperature and radius of the central 
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star. Constant Tq = '^^^^^ in unit of temperature, depends on stellar mass and radius R^,, 
and mass of hydrogen atom nip. 

Temperature and surface mass density for surface layer can be expressed as 

T. = (^)'t., (3) 

where a = ^ is the grazing angle of the irradiation from the central star to the surface of the 
flaring disk, and Hp{Ts) is the Planck mean opacity for surface layer. 

Assuming the majority of dust grains in the disk as silicates, with evaporation temperature 
Tevp — ISOOii', the disk will be truncated in a special radius Rrim, inside which, the grains are in 
the gaseous state mixed with gas of light elements and are optically thin to the stellar irradiation. 
In the case that self-irradiation cannot be ignored, the location and optical height of the rim can 
be jointly determined by following equations: 

\ ^ evp J \ i^rim / 
„ _ I T^evy \ 2 / Rrim \ ^ p 

\ Iq J \ R* J 

where the dimensionless constant Xrim is different from that of the interior layer of the flaring 
disk because of the different grain density distribution and disk structures. From the discussion in 
Garaud and Lin (2006), one can expect that Xrim will be some bigger than x, due to the optically 
thin inner hole. 

The above expression for Rrim is different from t he truncation radius for the disk gas. Man y 



T Tauri stars have strong (kG) fields on their surface (jJohns-Krull et al. 20041 : iJohns-Krull 20071 ) 



These fields dominate the gas flow and truncation protostellar disks out to a distance 

Rrn = /3//^M,-'/'(2GM,)-V7 (7) 

where fi is the dipole magnetic n ioment, is the accretion rate in the disk, and dimensionless 



parameter (3 ranges from 0.5 to 1 (jKonigl 19911 ). Since grains are the primarily opacity sources, the 
dust destruction radius is more relevant for radiation transfer and the location of Rrim provided it 
is larger than Rm- 

In most regions of the disk, especially at large radii, energy produced by the viscous dissipation 
is relatively small compared with stellar irradiation. For computational convenience, we consider 
the contribution of the frictional heat only in the shadowed region, stretching from the inner rim 
to the outer boundary Rji. The following result further justifies this simplification. By equating 
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with we obtain the locus of boundary at 



H 



R 



X 



2x1,, 



Rr 



(8) 



Following the framework developed by Ru den and Lin (jRuden and Lin 19861 ) on the basis of the 
a— disk model (jShakura and Sunyaev 19731 ) . we obtain the effective temperature Tg in the shadowed 
region as 



4amp J yAU'^J \AU 



where the dimensionless efficiency factor as — 0.001. 



3. X-Ray Ionization and Grain Charging 

Absorbed by grains and gas in the disk, the stellar X-ray can only penetrate a certain distance 
into the disk just like the visual photons do. Except for the directly exposed inner rim and the 
flaring surface layer, the structure of other regions in the disk is not significantly influenced by the 
X-ray irradiation. 



3.1. Ionization state of the gas 



The rate of change in ionization fraction of the gas is 



dx 
dt 



(l-x) 



r 



+ e(T)nx 



aB(T)nx H SeVeCre(Q)x 



(10) 



where T is the ionizing energy input rate in J s ^ H-atom ^, is the av eraged source temperature, 
asiT) = 3x 10~^^T~^/^ photons m^ s~^ is the case B recombination ra te (Osterbrock and Ferland 20061 ) 



e{T) = 6 X io-i7ji/2g-So/fcsT g-1 -g colhsional ionization rate (jPrain 2011 



and n 



IS 



the total number density of hydrogen. In the above equation, the last term takes into account the 
recombination of protons and electrons after they collide with the grains. Se, Ve, and de are the 
absorption rate, velocity dispersion, and grains' cross section for electron capture which may be a 
function of the grains' net charge (see subsections below). 

In the above equation, only photons that are more energetic than the ground state energy of 
the hydroge n atom. En = 13.6eV contribute to the flux of ioniz ing photons. The observed energy 
distribution (jHudson 199ll : ICrosbv et al. 19931 : IWolk et al. 20051 ) in the X-ray range suggests that 
T Tauri stars emit photon flux (ie photons per per second) 



dN^jE) _ No 
dE^ Eq 



E 

Eq 



-1.7 



(11) 
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where the normalization factor Nq can be obtained from the total X-ray luminosity Lx in the 
energy range between Eq and E^ax (which is 8kev for the COUP sources), 



Lx 



^^^^ 



iVoEb 
0.3 



Emax ^ 



En 



1 



(12) 



In equation (llOp . the heating and ionization input rate 



Eax{E) — ^T^dE 

Eo 



dE 



Ea 



(13) 



(14) 



where the X- ray photoionization cr oss section ax{E) = a^E/Eo) for photons' energy E, and 



a = IQ-^^m^ dGlassgold et al. 1997I '). From equations ([II|), ([l2|), and ([11]), we find 



r 



Npa 
3.2 



En 



En 



-3.2' 



7.5 X lO-^Lx? /E, 



RrimEo 



Ea 



-0.3 



(15) 



This expression can be used to determine the changing rate of ionization fraction in equation (jlOp . 



The internal thermal energy {Eg) equation of the gas can be written as (jLin and Murray 19921 ) 



rrir, 



dEg 

dt 



+ P 



d{l/p) 
dt 



T{l-x)-nA{x,T)-n 



-x^kBTasiT) + x{l - x)e{T)Eo 



(16) 



where p = pkT{\ + x)/mH is the pressure. At temperatures below 10^ K, the cooling function 
A(x,T) is given by Dalgarno & McCray (1972) for optically thin line emission by metal ions with 
a solar abundance. Similar to equation ()14p . we replace F with 



Nqo_ 
2.2 



1 



Er, 



En 



-2.2 



En 



En 



-0.3 



(17) 



In our model, the ionization fraction and the gas temperature of the gas are determined by the 
intensity of the ionizing photons emitted near the central star. During quiescence when the X-ray 
luminosity is less than ~ 10^'^ J s~^, the surface of the disk region at ~ 0.1 AU retains its low gas 
temperature and ionization fraction. In this state, contribution from the collisional ionization rate 
e(T) is negligible because free electrons and protons are more likely to collide with grains than 
directly with each other to recombine. With their relatively large cross section, grains easily can 
capture ionized electrons and protons, and become negatively charged. 

However, during some outbursts when the flux of X-ray increases to ~ 10^^ J s~^ Never- 
theless, Lx is generally small compared with the stellar luminosity at longer wavelengths such 
that Ti, Ts, and Tg are not affected by them. However, X-ray outbursts can significantly in- 
crease ionization fraction and the gas temperature, perhaps trigger the onset of thermal instability 
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(jMurrav and Lin 1992 ) in the tenuous upper layers the disk at ~ 0.1 AU. In this hmit, cohisions 
between the electrons and protons are important in the determination of the ionization fraction. 
Electrons' thermal energy is also elevated. 

Now let's start estimating the equilibrium temperature and ionization fraction of the disk gas 
during quiescent state and X-ray flares. We show that equilibrium can indeed be reached within a 
relatively short time compared with the daily timescale of X-ray flare duration. By taking the right 
hand side of both equation (jlOp and (jl6p to be zero, we can solve the equilibrium temperature and 
ionization fraction as Tg ~ 472Er and Xq ~ 0.038 in the quiescent state when the X ray luminosity 
Lx — 5 X 10^^ J/s. Note that the value of Tq is small compared with Tj produced by both viscous 
dissipation and reprocessed stellar photons in the visual wavelength range. Therefore the X-ray 
photons do not significantly modify the disk structure and ionization fraction during the quiescence. 

However, the equilibrium temperature and ionization fraction as Tj ~ 9700/^ and x/- ~ 0.79 
during the X-ray flares when the X ray luminosity increases to Lx — 1 x 10^^ J/s. In this phase, 
Tf > Ti and photoionization by the X-ray photons determines the thermal structure of the disk 
surface layer. In the above estimation, we are looking at gas/plasma with number density n ~ 
lO^^m"'^ at 6 scale heights from the midplane {xrim = 6) in the frontier of inner rim located 
around Q.lAU from the central star (refer to the numerical section). The optical radial thickness 
of the inner rim is around 0.04AC/ at 6 scale heights, which is comparable with Rnm ~ Q.lAU 
and Hrim ^ . Consistent with our previous intuition, both temperature and ionization get 

elevated during X-ray flares. 

For dynamics of gaseous ionization by X-ray flares, we flnd the ionization fraction goes as 
(l-e-*/^0 if we neglect the collisional ionization and direct recombination, with the response time 



r,= -V + ^^We ^lO^s (18) 



which is much shorter than one day. By taking into account the collisional ionization and direct 
recombination, the ionization should proceed even faster. For heating and cooling dynamics, we 
estimate the response time respectively from the equation (jl6p is of the order 

rn = MZI^.i, (19) 

rc = ^4^^-10. (20) 



n. 



KTf,Xf) 



which are again much shorter than the daily X-ray flare timescale. Finally we show that blackbody 
irradiation for grains is so effective that the conductive heat transfer from overheated ambient gas 
to grains is negligible. In fact, during X-ray flares, grains are heated not only by stellar flux, but 
also by the collisions with overheated gas. Thus instead of equation ([5]), we have 



'^T^rain = ^ ( 77^) ^ (l + |^ ) + Gnt-^A^B {Tf - Tgra^n) (21) 
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where a is the Stefan-Bohzmann constant, and average velocity of hydrogen atoms is approximated 
as VH = \J ^^^^ ■ For ambient gas/plasma temperature during X-ray flares Tj ~ 9700i^r, the equi- 
librium temperature for grains Tgyfii^i — 15007^!" which is comparable to the equilibrium temperature 
Tevp = 1500i(' in the quiescent state. Thus conductive heat transfer from ambient gas is entirely 
negligible compared with the stellar irradiation in determining temperature of grains in the inner 
rim. Last but not least, after the X-ray flare, the rise in the gas temperature implies that the disk 
scale height will increase by a factor of two on the local dynamical time scale (probably a few days) , 
which will further lower the local gas density. As a result, the equilibrium temperature should be 
higher, and more charge will be loaded on grains. Thus grains experience a stronger magnetic force 
and a smaller gas drag, but at the same time a outward gas push. It's complicated to judge which 
effect dominates, and we suspect all to be insignificant. So in this paper, we omit the influence of 
the change of gas scale height right after X-ray flares. 



3.2. Charge-loading process 

We now consider the charge carried by the grains. X-ray fl ares charge g rains mainly through 



two mechanism: collisional charging and photoelectric emission (iDrain 201 ll ). Photoelectric emis- 
sion happens when a X-ray photon excites an energetic electron out of the grain surface. However 
the effect of photoelectric emission is comparable with collisional charging only when the photo- 
electric yield is of order unity. Hindered by multiple complications, photoelectric yield is est imated 



around 0.1 for the grain size and photon energy we considered (jWeingartner et al. 20061 ). thus 



negligible. We omit the photoelectric emission and focus on the collisional charging effect below. 
In the framework of a cer tain grain with cha rge Q, the cross sections for electrons and protons 



to collide with this grain are (jUmebavashi 19831 ): 



2 ' 



(Je = 7rr^e*'^'=o -dBe (22) 



a„ = 7rr2(l--i (23) 

where e is the elementary electric charge, Eq is vacuum permittivity, Ef, = ^mefg and Ep = ^nipVp 
are the kinetic energy of protons and electrons respectively. As shown above, the temperature equi- 
librium can be achieved in a relatively short time, so that the energy should be equally partitioned 
between electrons and protons Ef, — Ep = ^ksTf. 

By introducing 5(e) and Sp as the absorption rate for the electrons and protons after their 
collision with grains. We write the microscopic charge balance equation for grains: 

SeHeVeCFe = SpUpVpap, (24) 

where n^, and Up are number density of protons and electrons respectively. By defining dimensionless 
parameter 

1 (25) 



47reo VdEe ' 
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we can rewrite ()24p as 

1 + e neSe\mp) ■ ^ ' 

High spee d electrons tend to h ave a lower absorption rate than that for protons, and we estimate 



0.01 ^ ^ < 1 (jUmebayashi 19831 ). Moreover as inferred from ionization fraction ~ 0.79, during 



X-ray flares, the electrons and protons are so abundant that the charge on grains consist of only 
a tiny portion of the entire plasma, thus we have ng = Up. With these restrictions, we find from 
equation (j26p 2.5 ^ < 6.4, which is entirely determined by the ratio of absorption rates for 
electrons and protons. We adopt 9 = 5 (with Sg = 0.05) in the following analysis. Conclusively we 
know in the equilibrium state, grains are negatively charged to 

Qo = -5 (27) 

e 

which is proportional to the radius of the grain and the electrons' energy E^, but independent 
of the ionization fraction and electron density. During X-ray flares, plasma temperature in the 
inner rim increases dramatically, leading to a increase in E^ from ^ksTq to ^ksTf, and eventually 
a prepositional increase in Qo by ^r- ~ 20 times. 

Dynamics of charge loading processes can be determined by the following differential equation 

^ = -eSeVeCFe{Q)ne + eSpVpap{Q)np (28) 
We can determine the response time for charge loading on grains to be 

TQ ~ , ° - lO^'s (29) 

which is very swift, thus it takes nearly no time for the grains to get charged. 

In conclusion, we show in this section that X-ray flares will not only significantly ionize gas 
in the inner rim of the circumstellar disk, but also elevate its temperature by thousands of Kelvin 
higher. Both ionization and temperature elevation will accomplish in a much short timescale 
compared with X-ray flare duration. The overheated gas and plasma are shown to have negligible 
effect on equilibrium temperature of the grains, but will significantly increase charge on the grains 
in a relatively short time. In next section, we will show that the newly charged grains will undergo 
a collective migration driven by stellar magnetic field, thus modulate the height of inner rim of the 
disk. 



4. Grain Dynamics 

When the temperature in the inner rim is high enough (~ lO^K) for coUisional ionization of 
metal atoms, magnetic field can be coupled effectively with the gas in the inner region of a circum- 
stellar accretion disk. However despite of voluminous studies on the interaction of stellar magnetic 
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field and the circumstellar disk (Pringle and Rees 1972[ : Ghosh and Lamb 1979a : 1979b : Alv 1980 : 
Konigl 199ll : iLivio and Pringle 1993 : IShu et al. 1994 : iHavashi et al. 19961 ). there is little consen- 
sus on how well the plasma couple with the field lines. Some outstanding issues include "whether 
the disk excludes the stellar magnetic field by diamagnetic currents or the field can penetrate and 
thread a large fraction of the disk, whether the threade d field rein ains closed (connecting the star 
and the disk) or becomes open by differential shearing (ILai 1999I') ." There has also bee n a recent 
trend on 3D-modeling magnetic stellar field (jGregory et al. 20071 : iJardine et al. 20071 ). which is 
very complicated and beyond the scope of our analytical framework. In this section, we first de- 
velop a parametric framework for the magnetic field, taking into account the interactions between 
the stellar field and the circumstellar disk; and then discuss the charged grain dynamics driven by 
the magnetic field. 

As mentioned above, the circumstellar disk would be disrupted at the magenetospheric bound- 
ary where the magnetic stress is large enough to remove the excess angular momentum of the 
nearly Keplerian fiow over a narrow transition zone and channel the plasma onto the stellar polar 
caps. The dynamic consequence of disk accretion onto the central star is to synchronize the stellar 
spin with the disk rotation. How much time it takes to reach an equilibrium state depends on the 
coupling int ensity of field line and the disk. The corotational radius Rco ~ Rm, and Konigl adopted 
Rco = '^Rm (jKonigl 199ll ). and in fact the ratio can be even smaller. Within the range from Rco to 
Rm the disk angular velocity departs significantly from the Keple rian value, and corotates with th e 
central star; while beyond Rco, the angular velocity is Keplerian (jChosh and Lamb 1979al : Il979bl ). 
For a typical T Tauri stars, the spin period is around 8 days, so that Rco ~ 0.07 AU, which as shown 
later, is smaller than Rrim the sublimation frontier of the disk. In our parametric framework, we 
denote 



Rc 



eRr 



(30) 



where the dimensionless parameter e ranges from 0.3 to less than 1. In order for the stellar magnetic 
field to induce a spin equilibrium (in which there is little net angular momentum flux between the 
stellar spin and the disk), it must be well coupled to the disk beyond the corotation radius (Ju et 
al in preparation). 

We assume the stellar field can influence the surface layer of the disk near Rrim- We consider 
the case that stellar magnetic dipole, stellar spin and the disk angular momentum are aligned with 
each other. In the disk mid-plane, the magnetic fleld in the cylindrical coordinate can be expressed 
as: 

i?2 



Br 



-B 



Rl 



0, 



(31) 

(32) 
(33) 



where Bg ~ lO^Gs represents the magnetic strength at the equator at the surface of the central 
star. The negative sign for the longitudinal component Bz comes from our consideration that 
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the magnetic moment is parallel to the stellar spin. We leave the alternative anti-parallel case 
in the below discussion section. The azimuth component B^p comes from twisting of B^, and 
the upper (lower) sign corresponds to the value at the upper (lower) disk surface. The quan- 
tity C specifies the azimuthal pitch of the field line. In general we expect C ^ 1; but its actual 
value or form depends on details of the dissipative processes involved in the disk-magnetic field 
interactions. I f the stellar magnetic f ield threads the disk in a closed configuration, we expect 



C oc {ujs — ^d) (ILivio and Pringle 19921 ). where ujg and ujd are angular momentum for central star 
and disk respectively . However, it has been argued that the differential shearing and the plasma 
flowing from the disk into the overlying magnetosphere will blow the field lines open and main- 
tain theni in an open configur ation, in which case we expect C to be positive and of order unity 



(jLai 19991 : iHavashi et al. 19961 ). The radial component Bji comes from the inward accretion flow 
with velocity vji ~ v/R which is much smaller than the azimuthal disk velocity, so we adopt = 
in the analysis. Because p''"" ^ 1 as shown in the following numerical example, when we discuss 
the grain dynamics, we can always approximate magnetic field as 

B = BdTCBs^^ (34) 
in the proximity of disk mid-plane, where B/j is the magnetic field by unperturbed stellar dipole. 

After grains in the inner rim get charged by X-ray ionization, they move upward (downward in 
antiparallel case) collectively driven by the magnetic field. As a consequence, both the shape and 
the location of the inner rim changes over time, with the scale of movement as the Lorentz radius 
Rl 

~ B^Qo' timescale as the Lorentz period Tl — We focus on the dynamics of 

grains which originally settle around the upper/lower brim of the inner rim, and try to infer the 
variability of Hrim based on these grains' migration. 

The equation of motion for a certain grain in the rotating framework of central star can be 
formulated as 

GM R, 

md^r = -rrid — ^;3^^ h QoVr X B - i/|v - vkKv - vk) - 2mdUJs x Vr - rudUJs X (^s X r), (35) 

where iVg is the spin angular velocity for central star, Vg = x R is the corotation velocity with 
the central star, Vr = v — Vg is grain's velocity in the rotating framework, and vk is the Keplerian 
velocity. The last two terms are Coriolis and centrifugal force in the rotating framework. We only 
preserve the radial component of gravity, because the vertical component has been balanced with 
the gas turbulence, and that's the reason why the grain can stay around the brim initially. We 
neglect collisions between grains. Gas drag coefficient u = ^vrr^p is calculated based on the fact 
that the velocity of the grain is much larger than the thermal velocity of the gas, so the gas can 



be treated as still relative to the fast grains. Hermite scheme (jAarseth 20031 ) can be adopted to 
numerically solve the differential equation, with constant time steps much smaller compared with 

_ 7 

the Lorentz period. Notice that the Lorentz force decreases as i? 2 ^ while the gravity decreases as 
R^^ . For the grains in the surface layer of the flaring disk around R > 3AU, the relative strength 
of the Lorentz force with respect to gravity becomes ~ 10~^ of the value at the inner rim. Thus 
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we can omit the magnetically driven dynamics of grains in the surface layer of the flaring region. 
After we get the grain dynamics, we approximate the variability of Hrim as the location variations 
of grains in the brim. 



5. Spectrum Variability 

The collective outward migration of grains (in the parallel case) puffs up inner rim of the disk. 
Consequentially, the irradiation flux from the inner rim increases; while larger region in the flaring 
disk becomes shadowed by the rim. We track the flux variability of both inner rim and shadowed 
region, by tracking the variability of the height of inner rim Hrim- As shown below, the blackbody 
irradiation for individual grains is so effective that it takes almost no time for the newly shadowed 
(exposed) grains to reach a new equilibrium temperature. So we can determine the boundary Rji 
dividing the shadowed region and the exposed region by enabling geometric relationship instantly, 
the same as the static case dS]): 

H{Rfi) _ H{Rrim) , X 

Once the variability of disk struct ure determined, computa tion of the spectra before and after onset 



of X-ray flares is pretty standard (iDullemond et al. 200ll ) 



Finally we justify the temperature response time for newly exposed/shadowed grains in the 
flaring part. For newly exposed grains in the surface layer, the temperature dynamics follows 

^rjaTH^r - AnrjeaT^ = l^4p,c,^. (37) 

For grains around 3AU, it's not hard to figure out the temperature response time is around 10""^ day, 
which is much shorter than the spectrum variability with daily timescales as we are interested in. 
Similarly, the newly shadowed grains in the surface layer undergoes a rapid drop of temperature with 
timescales around Is, much shorter than the daily timescale. As the irradiation from the surface 
layer can only penetrate into a distance in the interior layer, the newly exposed surface layer 
will heat a thin top interior layer quickly by the irradiation, without changing the temperature 
of the inner interior layer much. The surface mass density of quickly heated top interior layer 
ASj = Oj/ = — ~ — , then the temperature dynamics follows 

AaaTti^f - 2i^,aTt = -(-^ + Zc,)^, (38) 
R flH dt 

where = ^Rgas is the gas specific heat per mol. We estimate the temperature response timescale 
for grains in the interior layer around 2iAU around 10^ days, which is much longer than the daily 
timescale we are interested in. In the above discussion, we omit the effect of the vertically convec- 
tive heat transference. Since it has the time scale comparable to the local dynamics, the effect of 
convective heat transference should works at at timescale much longer than several days. Conclu- 
sively, the newly exposed/shadowed grains in the surface layer change instantly; while grains in the 
interior layer change little in the daily timescale. 
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6. Numerical Studies 

In this section, we apply our model to an numerical example to generate predictions and 
compare with observat ions. We take the central star as a spherical blackbody with = 1.8Mq 
dMuzerolle et al. 2009ll. = % = 4000 K. The outer cutoff radius of the disk Rout = 



2.3 X 10 AU (IChiang and Goldreich 19971 ). The radius of grains = 0.5fim. The value of x takes 



different in different regions, where y = 4 for the flaring disk and Yrim, = 6 for the inner rim due to a 



higher gas density (jPullemond et al. 200ll : iDullemond and Monnier 2010 : iGaraud and Lin 20061 ). 



The stellar magnetic field at surface Bg = lOOOGs, the azimuthal component coefficient in equation 
(f32]) ( = 0.5 and dimensionless coefficient in equation (f30l) e = 0.5. 

We first describe the whole physical processes of spectral variability triggered by X-ray flares. 
We assume stellar X-ray flare with luminosity ~ lO^^J/s starts at time t = 0, and lasts for one 
day. During X-ray outburst, gas in the inner rim gets ionized quickly, with the ionization fraction 
rising from around 4% to around 80% within 10^ seconds. Meanwhile the X-ray photons lead to 
a temperture enhancement of the gas in the disk surface layer from around AOOK to around lO^K 
within the same timescale. As a result, speed of electrons become higher, and the collisions between 
electrons and grains become more frequently. In the equilibrium state, which can be reached in 
almost no time, the charge on grains increases by around 20 times. Grains strongly negatively 
charged will be driven by the magnetic field to undergo collective migration inward, in the case 
that the moment of magnetic dipole is parallel with the disk angular momentum. Consequently, 
the inner rim of the disk gets squeezed, and the shadow casted by the grains on the flaring disk 
attenuates. As a result, the spectrum at shorter wavelength increases while the spectrum at longer 
wavelength decreases, as sketched in Fig. [TJ In our numerical example, the height of inner rim 
decreases by 17.6% in the first day, and decreases by 22.7% cumulatively in seven days after 
occurrence X-ray flares, as shown in Fig. [21 In fact as discussed before, location of the boundary 
of the shadowed region depends on the ratio of height over radial distance of the inner rim, we 
show variability of the height-radius ratio in Fig. [3l In Fig. HI we plot the trajectory of charged 
grains projected in the midplane. We find newly charged grains still undergo circular orbit, with 
a minor inward shift. Finally in Fig. [5l we show the variability reproduced by our model. A clear 
"seesaw" feature is shown here: the spectrum at shorter wavelength decreases and spectrum at 
longer wavelength increases over time, with a pivot around 8fxm. Spectrum at A = 5/xm decreases 
by 14.9% and 18.7% respectively within one day and the whole week; while spectrum at A = 15/im 
increases by 7.4% and 9.2% respectively within one day and the whole week. The result is consistent 
with observations described in the first section. 



7. Summary and Discussion 



In this section, we briefly discuss the influence of the direction and strength of stellar magnetic 
dipole, as well as the grain size on variability of the spectrum first, and then conclude the paper. 
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Fig. 1. — A sketch illustrating the spectral variability resulting from disk structural variability 
during X-ray flares. 




Fig. 2. — Height of inner rim decreases over time, in which case X-ray flare occurs at time t = 
day and ends at time t = 1 day. 
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Fig. 3. — Ratio of height over radius of the inner rim decreases over time, in which case X-ray flare 
occurs at time t = day and ends at time t = 1 day. 



In our numerical example above, the stellar magnetic moment is parallel with the angular mo- 
mentum of the disk. However, in a general situation, the angle between t he stellar magnetic moment 



and angular momentum of the disk can range from to vr, as discussed by (jFlahertv and Muzerolle 2010 ) 
Here we simply show the result for anti-parallel case (with the angle as vr). Contrary to the parallel 
case, grains strongly negatively charged will be driven by the magnetic field to undergo collective 
migration outward, so the inner rim of the disk gets puffed up. The variability of the shadowed 
region can be complicated: on one hand, the puffed- up inner rim cast larger shadow on the flaring 
disk; on the other hand, some stellar light may penetrate the attenuated puffed-up inner rim. As 
a result, the spectrum at shorter wavelength increases while the spectrum at longer wavelength 
may experience decreasing first and increasing later. The variability of the height of inner rim is 
shown in Fig. [6l and the trajectory of charged grains projected in the midplane is shown in Fig. 
[71 Because the gas gets much dilute outward, charged grains basically spiral out of the system. 

To look at the influence of magnetic field strength on variability of the spectrum, we perform 
the same numerical simulation except for stellar magnetic field with Bg = 500Gs and Bg = 2000Gs 
respectively in the parallel case. We expect stronger magnetic field means stronger Lorentz force, 
resulting in more variability. In fact, for stellar magnetic field with Bg = SOOGs, the height of inner 
rim decreases by 15.5% in the first day, and decreases by 20.6% cumulatively in seven days after 
occurrence X-ray flares. Spectrum at A = 5/.fm decreases by 13.4% and 17.2% respectively within 
one day and the whole week; while spectrum at A = Ibfim increases by 6.5% and 8.5% respectively 
within one day and the whole week. The numerical result indicates that weaker magnetic field can 
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Fig. 4. — Trajectory of charged grains projected in the midplane (in the reference framework 
corotating with the central star), in which case X-ray flare occurs at time t = day and ends at 
time t = 1 day. 



- 18 - 




A (yum) 



10 15 20 30 

Fig. 5. — Log- log plot of multiple spectra. The solid black line represents the total spectrum of 
the star-disk system at time t = day; the dashed red line represents the total spectrum at time 
t = 1 day, one day after the occurrence of X-ray flare; and the dot-dashed blue line stands for the 
spectrum at time t = 7 day, six days after the cease of X-ray flare. A clear seesaw feature is shown 
here: the spectrum at shorter wavelength decreases and spectrum at longer wavelength increases 
over time, with a pivot around Si^m. The spectral variability is around 20% within one week. 
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t(Days) 

7 

Fig. 6. — In the anti-parallel case, height of inner rim increases over time, in which case X-ray flare 
occurs at time t = day and ends at time t = 1 day. 



generate significant variability as well. On contrary for stellar magnetic field with Bg = 2000Gs, 
the height of inner rim decreases by 19.7% in the first day, and decreases by 24.9% cumulatively 
in seven days after occurrence X-ray flares. Spectrum at A = S/xm decreases by 16.4% and 20.1% 
respectively within one day and the whole week; while spectrum at A = Ibfim increases by 8.2% 
and 9.9% respectively within one day and the whole week. Spectrum indeed becomes more volatile 
for stronger stellar magnetic fleld. 

To look at the influence of grain size on variability of the spectrum, we perform the same 
numerical simulation except that grain radius = 0.1/im and = Ifim respectively in the 
parallel case. Mass of grains is proportional to r^, while the Lorentz force is proportional to Qq, 
which is proportional to r^, thus smaller grains experience larger force per unit mass, therefore are 
expected to display more variability. In fact, for grains with radius = O.l/xm, the height of inner 
rim decreases by 22.5% in the first day, and decreases by 21.7% cumulatively in seven days after 
occurrence X-ray flares. Spectrum at A = 5/xm decreases by 21.7% and 26.7% respectively within 
one day and the whole week; while spectrum at A = 15/im remains almost unchanged with decrease 
by 0.5% and 1% respectively within one day and the whole week. The numerical result is consistent 
with our analysis, and the spectrum also decreases at longer wavelength due to a sharp decrease 
in the spectrum from the inner rim. On contrary for grains with radius = 1/zm, the height of 
inner rim decreases by 15.5% in the first day, and decreases by 20.8% cumulatively in seven days 
after occurrence X-ray fiares. Spectrum at A = 5/im decreases by 11.9% and 14.8% respectively 
within one day and the whole week; while spectrum at A = Ibfxm increases by 8.1% and 10.7% 




Fig. 7. — Trajectory of charged grains projected in the midplane in the anti-parallel case (in the 
reference framework corotating with the central star) , in which case X-ray flare occurs at time t = 
day and ends at time t = 1 day. 
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respectively within one day and the whole week. 

In conclusion, we construct a physical model to attribute the puzzling rapid MIR spectral 
variability to the observed X-ray flares from T Tauri stars. Our simulation result captures all main 
observed features consistently, and is robust with respect to different grain size and strength of 
stellar magnetic field. We made various assumptions in the construction of this model, including 
the location of the Rrim and an assumed distribution of the disk aspect ratio. The latter can 
be mapped out with regular monitor of disk's SED changes. Both of these quantities may also 
depend on the stellar magnetic field. It may be useful to explore correlations, if any, between these 
observable quantities. Finally, the basic assumption that the SED variations are triggered by X-ray 
flares may be tested by simultaneous X-ray and infrared observations. 
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